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Extending chip performance beyond current 
limits of miniaturisation requires new materials 
and functionalities that integrate well with the 
silicon platform. Germanium fits these require¬ 
ments and has been proposed as a high-mobility 
channel material, 1 a light emitting medium in 
silicon-integrated laserspEl and a plasmonic con¬ 
ductor for bio-sensingP 5 Common to these di¬ 
verse applications is the need for homogeneous, 
high electron densities in three-dimensions (3D). 
Here we use a bottom-up approach to demon¬ 
strate the 3D assembly of atomically sharp dop¬ 
ing profiles in germanium by a repeated stack¬ 
ing of two-dimensional (2D) high-density phos¬ 
phorus layers. This produces high-density (10 19 
to 10 2 ° cm -3 ) low-resistivity (1CU 4 H em) metal¬ 
lic germanium of precisely defined thickness, be¬ 
yond the capabilities of diffusion-based doping 
technologies.® We demonstrate that free electrons 
from distinct 2D dopant layers coalesce into a ho¬ 
mogeneous 3D conductor using anisotropic quan¬ 
tum interference measurements, atom probe to¬ 
mography, and density functional theory. 

Doping germanium with homogeneous, free-electron 
concentrations well above the metal-insulator transition 
(ti 3 D = 10 17 cm -3 ) enables low-resistivity source/drain 
extensions in high mobility transistors,^ Ge-on-Si inte¬ 
grated lasers with maximal optical gainp^ and plasma 
wavelengths (2-5 um) suitable for biological sensing.® 
Mainstream top-down implantation is inadequate for 
this purpose because vacancy-enhanced diffusion and 
the formation of neutral complexes create electri¬ 
cal deactivation and doping profiles that are broad 
and inhomogeneous. 910 Self-limiting surface reactions 
provide a promising alternative to control doping 
processes from the bottom-up. 11 This approach has 
produced monolayer-doped semiconductors with high- 
density, strongly confined, two-dimensional electron 
gases (2DEGs)p2 Combinations of bot tom-u p and top- 
down approaches have been propose d 13 ^ to extend 
monolayer doping from 2D to 3D. While dopants are 


deposited in a single 2D layer, their distribution in 3D 
was obtained by thermal diffusion, with associated loss 
of atomic precision and profile homogeneity. 

Here we demonstrate an exclusively bottom-up ap¬ 
proach to produce an effectively 3D doped, high-density, 
low-resistivity metallic germanium with a precisely de¬ 
fined thickness and doping profile. This is achieved by 
the repeated deposition of N nearly-identical phospho¬ 
rus doped layers as illustrated in Fig. la. This approach 
preserves the vertical atomic-precision associated with 
monolayer doping and creates a homogeneous 3D sys¬ 
tem, provided the interlayer spacing is sufficiently small 
that the electrons can readily move between the layers. 
Each layer is prepared in a three-step process (Fig. la): 
(1) self-saturating chemisorption of phosphine molecules 
(PH 3 ) onto a clean Ge surface^® (2) substitutional in¬ 
corporation of P dopants into the Ge lattice by thermal 
annealing to provide a two-dimensional electron density 
n 2 _Dp (3) encapsulation with Ge by molecular beam epi¬ 
taxy to separate layers at a distance d. By stacking a 
large number of layers a precisely-doped slab of thickness 
h = (N—l)d is produced. Using an inter-layer separation 
d comparable to the Bohr radius of phosphorus in germa¬ 
nium («8 nrn), vertical electron delocalisation creates a 
homogeneous 3D system with a density = Nri 2 D/h. 
With current monolayer doping techniques 16 achieving 
ri 2 D hi the range of 10 13 to 10 14 cm -2 and with sub- 
10 nm layer separations, we can expect high 3D densities 
n^D in the range of 10 19 to 10 2 ° cm -3 . The transition 
from single 2D layers to an effectively 3D doped material 
is explored in this paper using a single-layer (Fig. Id), 
a bi-layer (Fig. le), and a multi-layer sample (Fig. If) 
grown in ultra-high vacuum using the sequence described 
in Fig. la (Methods section). We use an interlayer sepa¬ 
ration of d ~ 5.7 nm, which is less than the Bohr radius. 

Three-dimensional atom probe tomography of these 
samples (Fig. lc-e) shows dopant layers that are nar¬ 
row and well separated. The dopant distribution of the 
single-layer sample has a full width at half maximum of 
1.41 ± 0.05 nm. The P atoms are distributed randomly 
within the doping plane with a very high planar density of 
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FIG. 1. Three-dimensional assembly of atomically sharp doping profiles using a bottom-up approach, (a) 

Phosphorus doped layers in germanium are fabricated in ultra-high vacuum by adsorption of phosphine molecules (PHa) onto a 
clean Ge(001) surface, thermal incorporation of P atoms, and encapsulation under an epitaxial layer of germanium of thickness 
d. (b) Repetition of the sequence in (a) produces a highly doped Ge film of total thickness (N—l)d. Pulsed laser atom 
probe tomography results from a (c) single layer, (d) bi-layer, and (e) multi-layer (18 layers) samples showing the cross-section 
distribution of dopant atoms. 


1.44 x 10 14 cm -2 . A similar average width (1.4±0.1 nm) 
and density (1.2 ± 0.3 x 10 14 cm -2 ) is found for layers 
of the bi-layer sample. In the multi-layer sample the av¬ 
erage width and density per layer are 2.0 ± 0.4 nm and 
(1.2 ± 0.3) x 10 14 cm -2 , respectively. There is a gradual 
increase in width by 0.06 nm from one layer to the next 
due to the accumulated thermal budget of the repeated 
deposition process (see Supplementary Section 1 for de¬ 
tails). Crucially, we find that the inter-layer separation 
is preserved from layer to layer, averaging 5.65 nm with 
a variance of less than 0.05 nm. These metrics confirm 
that the vertical atomic-precision associated with mono- 
layer doping is largely preserved when multiple layers are 
assembled into a larger stack. 

By comparing atom probe tomography and Hall mea¬ 
surements (see Methods) we determine that approxi¬ 
mately 26 to 44% of the dopants are electrically active. 
Despite incomplete activation, presumably due to the for¬ 
mation of P-P dimers,^ the measured electronic planar 
densities Nn 2 D are high at 6.3, 12, and 56 x 10 13 cm -2 in 
the single layer, bi-layer and multi-layer (18-layer) sam¬ 
ple, respectively. This corresponds to ultra-high 3D elec¬ 
tronic densities of 4.5, 1.9, and 0.6 x 10 20 cm -3 . We 
measure exceptionally low resistivities of 2.0, 4.5, and 
6.7 x 10 -4 fl-cm that are consistent with the resistiv¬ 
ity vs. density dependence expected in heavily-doped 


bulk GeF 17 18 This demonstrates that our bottom-up, 
interface-free doping technology achieves bulk-like resis¬ 
tivities in extremely thin («1.4 nm) doping profiles. 

Quantum interference measurements at cryogenic tem¬ 
peratures and in a vector magnetic field ( B) allow us to 
probe electron motion both within and between dopant 
layers, and observe the evolution from 2D towards a ho¬ 
mogeneous 3D electronic system as the number of layers 
is increased. Figure 2 shows polar plots Act wl (i9) of the 
weak localization (WL) positive corrections to the mag¬ 
netoconductivity measured at constant field and variable 
angle ■&. The single layer sample (Fig. 2a) exhibits a 
strong anisotropy with A(t m L ('d) at a maximum when 
B is perpendicular to the dopant plane. As the mag¬ 
netic field is rotated, Aa WL ('d) collapses and is nearly 
negligible when B is parallel to the 2D layer. This find¬ 
ing is intuitive: the 2DEG is strongly confined in the 
vertical direction and hence the backscattered particle 
trajectories that give rise to WL are located within the 
2D dopant plane. In the bi-layer sample (Fig. 2c), the 
anisotropy is still present but diminished. This indicates 
that in addition to motion within the layers, electrons are 
able to jump between them and form coherence interfer¬ 
ence loops that are affected by the parallel magnetic field. 
The two layers are strongly coupled despite the layer sep¬ 
aration and effectively behave as a coherent 2D system 
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FIG. 2. Evolution of the electronic system from 2D to 3D probed by electronic transport in different magnetic 
field orientations. (a), (c), (e) Polar plots of the size of the weak localization positive corrections to the magnetoconductivity 
Act" L ( r i9) = ct xx ([ B| = 1 T, $) — cr xx (\B\ =0 T) as a function of angle i? at a temperature of 200 mK for the single layer, 
bi-layer, and multi-layer samples, respectively; the angle i? is defined with respect to the dopant plane, with the components 
of magnetic field perpendicular and parallel to the plane of the dopant layers given by B± = \B\ sin(d) and By = jB| cos($). 
Dashed grey lines are contours of constant Act" l ($). Black lines are theoretical fits to the data; (b), (d), (f) Temperature 
dependence of the zero magnetic held resistivity p xx (normalized to its value p xx at T = 5 K) measured for the same samples 
(lines are a guide for the eye), showing how the characteristic logT behaviour for 2D systems weakens as the number of doping 
layers is increased. Illustrations of the Hall bars and weak localisation backscattered particle trajectories are also shown for the 
three samples. 


of finite thickness. Quantitative analysis of the weak lo¬ 
calisation yields the time scales associated with elastic 
scattering (r e ), interlayer tunneling (r t ), and clephasing 
( 't v ). The analysis shows r t < r e <C t v (see Supplemen¬ 
tary Section 2 for details), confirming coherent tunnelling 
of electrons between layers over a timescale comparable 
to scattering off dopants within each layer. In the multi¬ 
layer sample, the polar plot (Fig. 2e) is nearly circular. 
The electron self-intersecting scattering paths are equally 
probable in all the three directions, supporting the pres¬ 
ence of a homogeneous 3D system. 

Temperature-dependent resistivity measurements at 
zero magnetic field provide additional confirmation of 
this 2D to 3D cross-over. Figures 2b, d, and f show how 


the logT dependence characteristic of a 2D system 1 ^ is 
increasingly suppressed as layers are vertically stacked. 
Overall, this suppression of the temperature dependence 
and the progressive loss of anisotropy of the electron 
quantum interference reflect the loss of vertical confine¬ 
ment as electron motion approaches that of a 3D metallic 
conductor. Despite the atom probe imaging demonstrat¬ 
ing that the sample comprises well-defined 2D layers, the 
strong inter-layer coupling means that the electrons see 
the sample as a bulk 3D material. 

The thickness h of the conducting 2D systems in the 
single and bi-layer samples can be extracted by fitting 
Act wl [$) in Fig. 2a, c to a gener alized angle-dependent 
Hikami-Larkin-Nagaoka expressio d 20 l 2 d (see Supplemen- 
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FIG. 3. Density functional calculations of stacked dopant layers in germanium, (a) Band structure of a single 
phosphorus dopant layer (6.3 x 10 13 cm -2 ) plotted from the 2DEG band minimum towards the zone centre. The bulk conduction 
band minimum (CBM) is indicated using grey shading, (b) Dopant potential and donor electron density for the same single 
layer structure, (c) Correlation of band energies for a pair of phosphorus dopant layers (both 6.3 x 10 13 cm -2 ) as a function of 
separation. The Fermi energy is indicated by empty circles. For reference, the separation in the experimental bi-layer sample 
(5.7 nm) is indicated by a vertical dashed line, (d) The dopant potential and donor electron density for the same pair of dopant 
layers at the experimental separation, (e) Dopant potential and donor electron density for an infinite stack of dopant layers 
(each 3.1 x 10 13 cm -2 ) with a spacing of 5.7 nm. 


tary Section 3 for more details) that includes an ad¬ 
ditional dephasing rate due to the parallel magnetic 
field . 22 2 31 We obtain thicknesses of h = 1.49 ± 0.03 and 
6.17 ± 0.05 nm for the single and bi-layer sample, respec¬ 
tively, in agreement with thicknesses of 1.41 ± 0.05 and 
7.3 ± 0.2 nm obtained by atom probe tomography. This 
confirms the interpretation of the multi-layered doped re¬ 
gion as a 3D space where electrons coming from 2D lay¬ 
ered dopants can freely move, thus realizing a 3D metallic 
conductor. 

The critical role of the interlayer spacing d on the elec¬ 
tronic transition from 2D to 3D behaviour is further ex¬ 
plored using density functional theory in which the acti¬ 
vated dopant densities are used as input parameters. A 
single phosphorus layer in germanium 24 is characterized 
by a pair of valley-split bands (Fig. 3a, labelled 11/ and 
21/), and a 2DEG density that is spatially confined to a 
width of « 7.3 nm by the self-consistent doping potential 
(Fig. 3b). Figure 3c describes the evolution of the band 
minima of two dopant layers as the separation between 
the layers is reduced. At large separation, the two lay¬ 
ers are effectively independent and the single layer band 
energies are preserved. At closer separation the band en¬ 
ergies split due to coupling between the two layers. This 
commences at d ss 10 nm and is clearly established at 
5.7 nm (dashed vertical line), the actual separation in 
our bi-layer sample. This coupling is also evident in the 
calculated electron density and doping potential (Fig. 3d) 
where the single layer densities and potentials are seen to 
overlap. In an infinite stack of layers (Fig. 3e) the over¬ 
lap between the electron density (and the potential) of 
each layer is further enhanced, supporting the finding of 
strong 3D inter-layer coupling in our multi-layer sample. 

In conclusion, our bottom-up approach to doping in 


Ge is capable of producing high electron densities (10 19 
to 10 2 ° cm -3 ) and low-resistivity ( 10 -4 fl-cm) metallic 
conductors of precisely defined thickness. As such, 
this technology has immediate relevance in electronics, 
photonics, and plasmonics, towards the development 
of high mobility transistors, industrially viable Si- 
integrated lasers, and rnid-IR plasmonics bio-sensors, 
respectively. Finally, tunable doping at high densities 
also provides an ideal test-bed to clarify open questions 
on doping-induced bottom-up superconductivity in 
group-IV semiconductors, as recently proposed by Shim 
and TahanPSl 


METHODS 

Sample preparation. All samples were fabricated in 
a customized ultra-high vacuum system (base pressure < 
5 x 10 -11 rnbar) comprising a MBE system (MBE Kom- 
ponenten) for Ge deposition and an additional chamber 
for surface preparation and PH 3 dosing. The chambers 
are connected via a UHV transfer tube. Ge(001) sam¬ 
ples 2.5 x 10 mm 2 in size were cleaved from a Sb doped 
Ge(001) 4 inch wafer (resistivity of 1-10 fl-cm). Atomi¬ 
cally flat, clean, and defect-free surfaces Ge(001) surfaces 
were prepared for all samples with the method detailed 
in Ref 1261 In brief, an ex-situ wet chemical treatment 
using HC1:H 2 0 (36:100) and H 2 02 :H 2 0 (7:100) is used 
to alternately strip and reform a germanium oxide pas¬ 
sivation layer. This is subsequently removed in-situ by 
a flash-anneal at 760 °C, followed by a 25-nm Ge buffer 
layer growth by MBE at a rate of 0.015 nm/s and sam¬ 
ple temperature of 400 °C. Prior to the first doping cycle 
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step, the surface is flattened by a final thermal anneal 
at 760 °C. For P doping, all samples were saturation- 
dosed at room temperature with PH 3 gas backfilling the 
UHV chamber at a pressure of « 5 x 10 -11 mbar via a 
leak valve. Thermal incorporation of P atoms into the 
Ge surface was obtained by increasing the sample tem¬ 
perature from room temperature to 400 °C at a rate of 
1 °C/s. Epitaxial growth of Ge spacers 5.7-nm-thick was 
performed by MBE with the following growth tempera¬ 
ture sequence. The first « 0.7 nm of the Ge spacer are 
deposited at 400 °C, the following ss 2 nm at 250 °C, and 
the final 3 nm, again, at 400 °C. Between the first and 
second step of the sequence, the growth is interrupted to 
allow for sample cool-down at a rate of 1 °C/s. This se¬ 
quence was engineered to minimize dopant diffusion and 
segregation whilst retaining a low-roughness surface at 
each doping cycle.®® The deposition process ends for all 
samples with a 30-nm-thick Ge cap layer obtained by 
extending the duration of the third step in the growth 
temperature sequence. 

Pulsed laser atom probe tomography. After re¬ 
moval from UHV, all samples were cleaved and one por¬ 
tion was used to produce atom probe compatible spec¬ 
imens. All samples were coated with protective amor¬ 
phous films of 25-60 nm Cr/20 nm Pt using a broad 
ion beam sputter system (Gatan 682 Precision Etching 
and Coating system). The coated samples were then 
transferred to a focused ion beam (FIB)/scanning elec¬ 
tron microscope (SEM) dual-beam system (FEI Nova 200 
Nanolab, Hillsboro, OR). Atom probe compatible needle 
specimens were prepared in the dual beam system us¬ 
ing site specific lift-out techniques, then mounted to Si 
microtip posts (CAMECA Atom Probe Technology Cen¬ 
ter, Madison, WI) and finally annular milled. Pulsed 
laser atom probe tomography was performed at the Uni¬ 
versity of North Texas Center for Advanced Research 
and Technology (CART) using a Local Electrode Atom 
Probe (LEAP) 3000X HR (CAMECA Atom Probe Tech¬ 
nology Center, Madison, WI) laser pulsed local electrode 
atom probe with a reflectron lens. Samples were anal¬ 
ysed at a base temperature of 30-50 K in laser pulsed 
field evaporation mode using a pulsed laser with a wave¬ 
length of 532 nm, pulse width of 10 ps, applied at a 
pulse frequency of 16 kHz, laser energy of 0.2-0.3 nJ 
and an evaporation rate of 0.001-0.005 ions/pulse. The 
resulting tomographic atom probe data was analysed us¬ 
ing the atom probe reconstruction software, IVAS 3.6.1 
(CAMECA Atom Probe Technology Center, Madison, 
WI). Reconstructions were correlated with Transmission 


Electron Microscopy and Secondary Ion Mass Spectrom¬ 
etry data (see Supplementary Section 1 for more details). 

Electrical characterization. Trench-isolated Hall 
bars structures to investigate the electrical properties of 
the doped layers were defined by a CHF 3 /CF 4 based dry 
etch with thermally evaporated A1 Ohmics connecting in 
parallel all multiple P layers of the doped stack. Electri¬ 
cal characterization at 4.2 K was performed using a dip¬ 
stick in liquid helium equipped with a superconducting 
magnet providing a perpendicular magnetic field up to 
2 T. Characterization of the device at lower temperatures 
between 0.2 and 5 K was performed in a cryogen-free 
dilution refrigerator, equipped with a triple axis vector 
magnet system (Leiden Cryogenics B.V.). This enables 
independent control of both perpendicular and parallel 
components of the magnetic field, with respect to the 
average dopant plane, and allows for magnetic field rota¬ 
tions to be performed at fixed field. The vector magnet 
was critical to performing reliable WL measurements to 
extract the 2D layer thickness since the large anisotropy 
in Act wl requires alignment of B to better than 0.5 de¬ 
grees with respect to the doping plane. For all measure¬ 
ments we used a four-probe setup using standard low 
frequency lock-in techniques and low injection currents 
(« 1 nA) to measure simultaneously the magnetic field 
dependence of the longitudinal p xx and transverse p xy 
components of the resistivity tensor, where x and y are, 
respectively, the directions parallel or perpendicular to 
the current flow in the Hall Bar. The longitudinal a xx 
and transverse a xy (Hall) conductivity were calculated 
from the measured resistivities via tensor inversion. 

Density functional theory calculations. DFT cal¬ 
culations on stacked phosphorus dopant layers in ger¬ 
manium were conducted using the SIESTA software 28 
and methods described for single Ge:P layers in Ref. [24j 
The DFT equations were solved using an atom-centered, 
double-numerical-plus-polarization (DNP) basis set and 
the local density approximation (LDA) with empirical 
on-site (+U) correction. The single and double dopant 
layer structures (Fig. 3a, b and Fig. 3c, d) were rep¬ 
resented using highly elongated germanium super cells 
of 300 atomic layers which is sufficient to separate the 
dopant layers from their periodic repeats. For the re¬ 
peated dopant layer stack (Fig. 3e) a much smaller unit 
cell of 40 layers was used in order to match the experi¬ 
mental layer separation of 5.7 nm. Phosphorus densities 
of 6.3 x 10 13 cnR 2 and 3.1 x 10 13 cm -2 in the dopant 
plane were represented using the mixed-atom approach 
described in Ref. [Ml All dopant atoms in the calcula¬ 
tions are confined to a single atomic plane. 
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SI. ANALYSIS OF PHOSPHORUS DEPTH PROFILES OBTAINED BY ATOM PROBE TOMOGRAPHY 

Atom probe tomography was used to create three-dimensional tomographic reconstructions of the single layer, bi¬ 
layer, and multi-layer samples which are shown in Fig. lc-e of the main text. In Figures Sla-c atom probe tomography 
data is integrated into phosphorus atomic concentration profiles and compared to concentration profiles obtained by 
Secondary Ion Mass Spectrometry (SIMS). This comparison highlights the superior depth resolution of atom probe 
tomography to characterise the sharp dopant profiles in our samples. The lower resolution in the SIMS data is due to 
ion beam mixing from the sputtering process, which is known to lead to artificial peak broadening. The atom probe 
phosphorus profile for the single layer sample in Fig. Sla is slightly skewed to the left (i.e. towards the sample surface), 
which indicates that a degree of dopant segregation during growth is present. We fit this profile with an exponentially 
modified Gaussian (black line in Fig. Sla), which describes a diffusion-broadened profile with an exponential leading 
edgeP 
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We use as fitting parameters the 1/e decay length A, and the amplitude A, position fi, and width w of the 
deconvolution Gaussian. We find A = 0.71 ±0.03 nm, which corresponds to a segregation length scale of approximately 
five atomic layer in the germanium crystal. This shows that the growth process induces very little dopant segregation 
and hence the doping profiles are sharp at the atomic level. 

The peaks in the bi-layer and multilayer concentration profiles (Fig. Sib, c) appear symmetric. The small inter¬ 
layer spacing limits phosphorus segregation in these samples and the peak shape is primarily determined by dopant 
diffusion from the accumulated thermal budget of the repeated deposition process. Hence, in these samples the doping 
profiles are better fitted with Gaussian functions (black lines in Fig. Sib and c). 


Single layer Bi-layer(average) Bi-layer (peak 2) Bi-layer (peak 1) Multi-layer (average) 
FWHM (nm) 1.41 ±0.05 1.4 ±0.1 1.47 ±0.05 1.32 ±0.03 2.0 ± 0.4 

Area (10 14 cm -2 ) 1.44 1.2 ±0.3 1.47 1 1.2 ±0.3 


Table SI. Analysis of phosphorus concentration profiles for single layer, bi-layer, and multilayer sample shown in Fig. Sla-c. 

Table SI summarises for all samples the key parameters reported in the main text. These parameters are the 
peak full width at half maximum (FWHM) and the area under the peak, obtained by numerical integration. In 
the single layer the FHWM is indicative of the dopant layer thickness h. We assume for FHWM in the single 
layer the same percentage error (3.6%) of the deconvoluted Gaussian FWHM obtained from fitting Eq. (1) to the 
data. In the bi-layer, the thickness is estimated by adding the finite average peak width to the interlayer separation, 
h = d ± FHWM = 7.3 ± 0.2 nm. For all samples, the area under each peak is indicative of the planar P density per 
layer (np). The values achieved are comparable to previous studies of single monolayer doped Ge:P sanrplesP 

In the multilayer sample, the average peak FWHM is increased to approximately 2 nm. To investigate the origin 
of this broadening, the FWHM for each layer are shown in Fig. Sle. The data shows that the width progressively 
decreases from the first deposited layer to the last. This can be understood in terms of an accumulated thermal budged. 
The first deposited layer is the broadest because it undergoes a small amount of thermal diffusion during the annealing 
and growth step of each subsequent doping cycle. The last deposited layer is the sharpest because experiences only a 
single thermal processing step. A linear fit of FWHM versus Doping Cycle quantifies this broadening as approximately 
0.06 nm per doping cycle. Despite this broadening, the layers are still distinct and their separation preserved. This 
is illustrated in Fig. Sle, where the linear increase in the peak position with the doping cycle yields an inter-layer 
separation value of 5.65 ± 0.05 nm, as reported in the main text. 


S2. QUANTITATIVE ANALYSIS OF THE WEAK LOCALISATION FEATURE. 

In this section we provide further details into our measurements of the magneto conductivity a xx (B) with orientation 
of the vector magnetic field perpendicular (B±) or parallel (£?m) with respect to the 2D doping layers. We perform 
a quantitative analysis of the weak localisation correction in perpendicular magnetic field to obtain the relevant 
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Figure SI: Analysis of atom probe tomography results. Phosphorus atomic concentration profiles from a (a) single 
layer, (b) bi-layer, and (c) multi-layer (18 layers) obtained by atom probe tomography (coloured circles) or SIMS (white 
circles) as a reference. Black lines are fit to the data. Analysis of individual peaks of the multi-layer sample showing: (d) 
planar doping density by integrating the area under each peak (dashed line is the average value); (e) full width at half 
maximum (FWHM) obtained by a Gaussian fit (dashed line is the average value);(f) peak position. Black lines in (e) and (f) 
are linear fits to the data from which we obtain the average FWHM broadening and inter-layer separation, respectively. In 
(d)-(f) the first doping cycle corresponds to the peak 1 in the concentration profile in (c). 


timescales associated with electron transport, indicating strong electronic coupling between stacked dopant layers. 
Throughout the range of temperatures investigated, electron transport is in the diffusive regime with very short 
measured elastic scattering times r e of 0.9 x 10~~ 14 , 1.0 x 10 -14 , and 1.8 x 10~ 14 s in the single layer, bi-layer, and 
multi-layer sample, respectively. 

Figure S2a, b, and c compare the conductivity a xx as a function of B± and Bn (black and grey curves, respectively) 
for the single layer, bi-layer, and multi-layer, respectively. In the single layer (Fig. S2a) weak localisation due to phase 
coherent backscattering reduces significantly the conductivity at zero magnetic field. The characteristic signature of 
weak localisation is a peaked positive magnetoconductance in perpendicular field, as B± threads flux through closed 
loops of particle backscattering trajectories and kills the weak localisation feature. In contrast, when the magnetic 
field is parallel to the 2D doping plane, very little flux is threaded through the particle trajectories, so almost no 
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suppression of weak localisation is observed, and, consequently, the magnetic field dependence is weak. At B = 1 Tesla 
(circles in Fig. S2a) a strong anisotropy is established. This anisotropy is also evident in the polar plot in Fig 2b of 
the main text. In line with previous weak localisation measurements on Ge:P single layersp we observe in parallel 
magnetic field a small negative weak anti-localization correction, due the local moments in the presence of strong 
Coulomb interactions. 

In the bi-layer sample (Fig. S2b), the perpendicular field dependence resembles that of the single layer sample. In 
the parallel direction, there is now significant field dependence because Bn threads flux through particle backscattering 
trajectories formed in-between layers due to strong coupling. In this sample, the anisotropy at 1 Tesla is still present 
but diminished. In the multi-layer sample (Fig. S2c), the perpendicular and parallel magnetoconductance are seen to 
be almost overlapping, and hence the WL correction to the magnetoconductivity is nearly isotropic. 

In Figure S2cl we show the perpendicular field magnetoconductivity per layer cr xx {Bj_)/N with TV = 1, 2, and 
18 in the single-layer, bi-layer and multilayer sample respectively. The measurements are over an extended range 
of magnetic fields to compare directly the magnitude of the WL correction across the three samples. We note a 
progressive suppression of the WL feature around zero magnetic field as the number of layers is increased, confirming 
the evolution of the electronic system from 2D to 3D. Note that if the layers were fully decoupled, all samples 
would show a similar dip in the magnetoconductance per layer. Associated with the dimensional cross-over is also 
the approximately twofold increase in electron mobility from the single-layer to the multi-layer. This is due to the 
population of 3D states extended in the vertical direction with reduced scattering from impurities within the doping 
planes. 

The relevant timescales of electron transport in the single layer and an insight into electronic coupling in the 
bi-layer sample are obtained by fitting of the WL feature in perpendicular field. For all fits we have subtracted 
the classical correction to the Drude conductivity a D {iiB) 2 . In the single layer (Fig. S2e), the magnetoconductivity 
Aer xx = cr xx (B±) — cr xx (0) is fitted to a phenomenological generalized Hikami-Larkin-Nagaoka expression for the 
quantum interference correction to the conductivity!^ 
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In this expression a and /3 are positive constants close to unity, F(x) = F(0.5 + x) — ln(x) with F{x) being the 
digamma function, D is the two-dimensional diffusion constant, t s and are the timescales characterising quasielastic 
spin scattering and dephasing, respectively. Excellent agreement between fit (black curve) and experimental data is 
obtained over the entire range of magnetic fields using values for the fitting parameters a, /3, r s , and t v listed in table 
S2. Note that fitting with a simple Hikarni model (Fig. S2e, dashed line), recovered by assuming /3 = 0 in Eq. (2), 
neglects quasi elastic spin-scattering and fails to describe the data over the entire range of the magnetic field. In the 
bi-layer sample (Fig. S2f), excellent agreement between fit (black curve) and experimental data is also obtained using 
Eq. (2), i.e. assuming that the two layers are coupled such that they effectively behave as a single coherent layer. To 
gain further insight into the electronic coupling in the bilayer, we also present in Fig. S2f theoretical fits (grey curves) 
to the expression: 
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This expression describes tunnel-coupled double quantum wells with a characteristic time r s as in Ref [J, and is 
generalized to include quasi elastic spin scattering as in Eq. (2). For strong tunnel coupling (r t < r e -C t v ), Eq. (3) 
reproduces the black curve in Fig S2f, recovering the single layer limit. When the layers are weakly coupled (grey 
curve in Fig. S2f; r 4 /r e = 10) or fully de-coupled layers (light grey curve in Fig. S2f; rt/r e = 1000) the theoretical 
fits fail to fully reproduce the experimental data. Overall, this quantitative analysis confirms that electron transport 
in the bi-layer is in the strong coupling regime. As described in the main text, this regime is characterised by the 
coherent tunnelling of electrons between layers over a timescale comparable to the scattering off dopants within each 
layer. 


S3. FITTING PROCEDURE TO EXTRACT LAYER THICKNESS FROM THE ANGULAR 
DEPENDENCE OF THE WEAK LOCALISATION CORRECTIONS TO THE 

MAGNETOCONDUCTIVITY 


The thickness h of the samples showing two-dimensional character, i.e. the single and bi-layer, is extracted by 
a quantitative analysis of the angular dependence of the weak localisation corrections to the magnetoconductivity 
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Figure S2: Magnetoconductivity in perpendicular and parallel magnetic field. Conductivity as a function of 
perpendicular or parallel magnetic field (black and grey lines, respectively) for a (a) single layer, (b) bi-layer, and (c) 
multi-layer (18 layers) at a temperature of 200 mK. (d) Conductivity per layer as a function of perpendicular magnetic field 
for the single layer, bi-layer, and multi-layer sample obtained by dividing the black curves in (a), (b), and (c) by 1, 2, and 18, 
respectively, (e) Weak localization correction to the magnetoconductivity for the single layer sample and theoretical fittings 
(black curves) obtained with the models described in the text, (f) Weak localization correction to the magnetoconductivity 
for the bi-layer sample and theoretical fits (grey and black curves) obtained for different electronic coupling regimes. Models 
are described in the text. 


A a WL (d) = <t(|.BJ, $) measured at \B\ = 1 Tesla. As reported for 2DEGs hosted in other materials systems, such 
as Si MOSFETsp Si:P delta doped layers,® and graphenep the effect of a parallel magnetic field is to enhance the 
electron dephasing rate by 


p* 


r -i + T -! 

'in l 1 Rn 


According to Ref. [5J the additional dephasing rate in parallel magnetic field is given by 


( 4 ) 



( 5 ) 


where l e is the mean free path (obtained from resistivity and Hall effect measurements), Z is the root mean square 
amplitude of the surface describing the two-dimensional electron gas, and R is correlation length over surface height 
fluctuations. In our angular measurements performed at unitary magnetic field B^ = sin 2 t?. As in RefUJ we assume 
that Z and R correspond to the thickness h of the doping distribution probed electrically and the mean donor spacing 
within a layer, respectively. R is estimated from the atom probe tomography analysis as 1 /yjnp. We extract h by 
fitting Aer WL (i?) (black lines in Fig. 2a-c) to the phenomenological generalized Hikami-Larkin-Nagaoka expression 
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Single layer 

Bi-layer 

a 

1.075 ±0.004 

1.056 ±0.007 

p 

0.617 ±0.001 

0.693 ±0.001 

T e 

0.0097 

0.01 

T V 

25.0 ± 0.2 

30.0 ± 0.04 

T s 

0.61 ±0.01 

0.34 ±0.01 

Tt 


< 0.01 


Table S2. Transport parameters for the single layer and bi-layer sample. r e was obtained directly by resistivity and Hall 
effect measurements. All other parameters were obtained by fitting the weak localisation feature in the magnetoconductance 
at a temperature of 200 mK. 


given in Eq. (2) with B± = cos'd for \B\ = 1 Tesla and using the effective dephasing time r^* defined in Eq. (4) 
instead of t v . Z is used as a single fitting parameter and is linked to the effective dephasing time As reported 
in the main text, we obtain thicknesses of 1.49 ± 0.03 nm and 6.17 ± 0.05 nm in the single and bi-layer sample, 
respectively. These values are both in agreement with the thickness obtained by atom probe tomography, confirming 
that the conducting regions in our samples match the extent of the stacked doping profiles. 
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